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We argue that the hghtest supersymmetric particles (LSP) can be pro- 
duced with extremely high energies E > 10^° GeV in the Universe at 
the present epoch. Their most probable sources are decaying superheavy 
particles produced by topological defects or as relic Big Bang particle. We 
discuss the mechanisms of production of LSP at ultra-high energies (UHE) 
and the interaction of the UHE LSP with matter. The most attention is 
given to the neutralino as LSP, although the gluino is also considered as 
a phenomenological possibility. 
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1 Introduction 



Cold Dark Matter (CDM) is probably the most abundant form of matter in the 
Universe. A natural CDM candidate is the lightest supersymmetric particle 
(LSP) which is stable, if i?-parity is conserved. In this Letter, we address the 
question if the LSP can be also among the ultra-high energy (UHE) particles 
filling the Universe. Theoretically the best motivated candidates for LSP are 
the neutralino and gravitino. We shall not consider the latter, because it is 
practically undetectable as UHE particle. Therefore, we pay most attention 
to the neutralino. 



In all elaborated SUSY models the gluino is not the LSP. Only, if the dimension- 
three SUSY breaking terms are set to zero by hand, gluino with mass rUg = 
0{1 GeV) can be the LSP There is some controversy if the low-mass win- 
dow 1 GeV < rrig < 4 GeV for the gluino is still allowed . Nevertheless, we 
shall study the production of high-energy gluinos and their interaction with 
matter being inspired by the recent suggestion (see also [^]), that the atmo- 
spheric showers observed at the highest energies can be produced by colourless 
hadrons containing gluinos. We shall refer to any of such hadron as ^-hadron 
(G). Light gluinos as UHE particles with energy E > 10^^ eV were considered 
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in some detail in the literature in connection with Cyg X-3 [^0. Additionally, 
we consider heavy gluinos with rrig > 150 GeV f^. 

Three mechanisms for the production of UHE LSP can be identified. 

(i) They can be produced in " astrophysical accelerators" due to the inter- 
action of accelerated protons with ambient gas. This mechanism effectively 
works only in the case of a light gluino, and it was exploited in 80s during 
the Cyg X-3 epic where the glueballino {gg bound state) was one of the main 
characters PJ^l- -^"^^ UHE ^-hadrons this production mechanism was indicated 
in Ref. 0]. The main difficulty of this mechanism is the low LSP flux, which is 
caused by the small cross-section for gluino production if rrig > 1 GeV and by 
the low density of the target nucleons or photons around the " accelerator" . 

(ii) Evaporating black holes are another possibility. High energy particles can 
be produced during the final stages of evaporation; the resulting spectra of 
cosmic rays have been discussed in Refs. @]. However, when these spectra 
are combined with various observational bounds on the mass fraction of the 
universe in black holes (see |0 for a review), one finds that the UHE CR flux 
from black holes is well below the observed flux. 

(iii) Decays of the supermassive particles produced either by topological de- 
fects or by Big Bang as relic particles, can naturally provide the large fluxes 
of UHE LSR 

The plan of this paper is as follows: We first discuss the possible sources of 
UHE LSP. Then we examine the cascade-production of LSP and their resulting 
spectrum. After that, we calculate the LSP fluxes for the two most promising 
sources. Finally, we examine the interactions of the UHE LSP with matter 
and discuss the status of (quasi-) stable gluino. 



2 Topological defects and supermassive relic particles 



Topological defects and supermassive relic particles are the two most promis- 



ing sources. Topological defects |T^ such as superconducting strings, monopoles, 
and monopoles connected by strings can produce UHE particles [|ll| . Here we 
will concentrate on cosmic necklaces, since this model seems to provide the 
largest UHE particle flux for fixed density of electromagnetic cascade radia- 
tion. 

Cosmic necklaces are hybrid defects consisting of monopoles connected by a 
string. These defects are produced by the symmetry breaking G ^ HxU{l) ^ 
H X Z2. In the first phase transition at scale rjm, monopoles are produced. At 
the second phase transition, at scale r]s < rjm, each monopole gets attached 
to two strings. The basic parameter for the evolution of necklaces is the ratio 
of monopole mass and the mass of the string between two monopoles, fid, 
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where /x ~ 77^ is the mass density of the string and d the distance between two 
monopoles. Strings loose their energy and can contract due to gravitational 
radiation. As a result, all monopoles annihilate in the end producing super- 
heavy Higgs, gauge bosons and their supersymmetric partners which we call 
collectively X-particles. The rate of X-particle production can be estimated 
as 

dnx r'^fJ^ . 

^ . i 

dt t^mx 

Similar to the case of UHE protons, the flux of UHE LSP is determined mainly 
by two parameters, r^/i and mx, which values must be of order 10^''' GeV^ and 
10^^ GeV, respectively, to have the flux close to the observed one. For a more 
complete discussion see Ref. [|12|. The diffuse flux of LSP produced by the 
decay of X-particles from necklaces is given by 

hsp{E) - - R{E) , (2) 

where dNisp/ dE is the spectrum of LSP from the decay of X-particles. 
Furthermore, R{E) is the attenuation length \{E) of the LSP if X{E) < do 
and R{E) = do otherwise, where to is the age of the Universe. 

If the neutralino x is the LSP, X{E) is much larger than do- It is determined 
as X{E) = {l/E ■ dialogs/ dt)~^ and the dominant contribution is given by the 
scattering of neutralino off background neutrinos, x+^bb ^ X + ^- The largest 
cross-section, a ~ G'jp{E/m^yel, where is the energy of the relic neutrino, 
is provided by the Higgsino component of the neutralino with Z'^-exchange 
in the t-channel. Apart from the smallness of the cross-section, A is further 
reduced by the small energy transfer in one collision. The Universe becomes 
transparent for UHE neutralino at red shift z < 1 ■ 10^, however production 
of neutralinos at early cosmological epochs is not important in the case of 
topological defects. 

The dominant energy-loss process of the ^-hadron is pion production in colli- 
sions with microwave photons. Pion production effectively starts at the same 
Lorentz-factor as in the case of the proton. This implies that the energy of 
the GZK cutoff is a factor rrig/mp higher than in case of the proton. The 
attenuation length also increases because the fraction of energy lost near the 
threshold of production is small, ^/mg, where /i is a pion mass. Therefore, 
even for light ^-hadrons, nig > 2 GeV, the steepening of the spectrum is less 
pronounced than for protons. 

Let us now come to superheavy relic particles as a source of UHE LSP. This 
case is formally identical to the production of X-particles by topological de- 
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fects if one replaces the RHS of Eq. (|l]) by nx/Tx- Following Ref. [|T3 



we 



will use the ratio rx = ^x^o/^x, = Px/pcbm as parameter characterizing 
the model (here px and pcDM are the mass density of X-particles and of the 
total CDM, respectively). The flux can be calculated by Eq. (|]), where now 
R{E) is the size of galactic halo Rh- We adopt the following astrophysical 
parameters for this case [|l^: Rh = 100 kpc, ^cdm = 0.2h'^ and h = 0.6. We 
assume for the halo density rnxUx = ^xPcdm and for the extragalactic density 
mxT^x = ^x^cdmPct, where ^x describes the fraction of X-particles in CDM 
and ficDM is the CDM density in units of the critical density pcr- 



3 Cascade-production of LSPs at the decay of X-particle 



The decay of X-particle results in a particle cascade similar to the QCD cas- 
cade in e~^e~ annihilation. The basis of the cascade development in both cases 
is given by probability p of production of an extra parton, p ^ g'^lnQ'^ > 1, 
where 5^ is a coupling constant. In the process of the cascade development, the 
energy and the virtuality of the cascade particles diminish progressively. 
Until the virtuality Q remains larger than the SUSY scale M|usy ~ (1 TeV)^ 
the decay channels to the usual particles and their supersymmetric partners 
have equal probabilities and the number of SUSY particles at each genera- 
tion is exactly equal to that of usual particles. When reaches Mg^gy, the 
supersymmetric particles go out of equilibrium and decay to the LSP. 

We performed a simplified Monte-Carlo simulation^ of the cascade includ- 
ing as elementary processes the transition probabilities between fermions, 
sfermions, gluons, gluinos, VT-bosons and winos, photons and neutralinos with 
probabilities similar to that in QCD. For simplicity, we assumed a common 
mass MsusY = 1 TeV for all SUSY particles except the LSP. We followed 
the evolution until a particle reaches the virtuality < Mgugy. In the case 
of the neutralino as LSP, each SUSY particle is then turned into neutralino 
after one or several decays. The case of heavy gluino, rrig > 100 GeV is very 
similar, with the obvious difference that gluino is turned into ^-hadron at the 
confinement radius. 

The case of a very light gluino is different. After SUSY particles go out of 
equilibrium, gluinos still participate in the QCD cascade due to g g + g. 
The processes of gluino production and radiation of gluons by gluinos are very 
similar to that of quarks and one can expect that in the low-energy part of 
the cascade the number of gluinos is roughly equal to the number of quarks. 
However, in the problem considered here, we are not interested in particles 
with too small x and our simulation does not include this low-energy regime. 



A more complete Monte-Carlo simulation of the cascade is in preparation |14|. 
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The obtained spectrum of the LSP can be well approximated for the energies 
at interest E ^ MgusY by 



dA^ I f E \-^-^ 



X 



with k ~ 0.25 for both neutralino and gluino as LSP. The fraction of energy 
transferred to the LSP is /susy ~ 0.4. In Fig. 1 and 2 we show the result- 
ing LSP fluxes from decaying X-particles and cosmic necklaces, respectively. 
Moreover, the neutrino fluxes calculated in MLLA approximation and exper- 
imental data |]16| are shown. 



4 Interaction of UHE neutralino with matter 



Let us consider the interactions of the neutralino x relevant for their detection. 



They are somewhat similar to the calculations ||T5[ for a photino. Mainly two 
processes are important for its interaction with matter, namely the neutralino- 
nucleon scattering x + — >^ all and resonant production of selectron off elec- 
trons X + e ^ e all. The first process is based on the resonant subprocess 
X + q ^ q sl\ and on neutralino-gluon scattering. The latter subprocess 
is important, because for high energies, and consequently for small scaling 
variable x, the gluon content of the nucleon increases fast. 

The resonant cross-section of the parton process x ~^ q ~^ qL.R all is given 
by the Breit-Wigner formula as 

Pern 1-5 ^^■^L,r) + ^^^L,i?J- tot 



where s = 2E^mN, E^ denotes the energy of the incident neutralino, uin 
the nucleon mass, Ml^r are the masses of left- and right-chiral squarks and 
s = sx. The total decay width Ftot of the squark depends strongly on the 
mass spectrum of the model. In the following, we parameterize our ignorance 
by Ptot = zV{qi R ^ q + x) with z > 1. To obtain the neutralino-nucleon 
cross-section, the parton cross-section a^^^{s) has to be integrated with quark 
distribution functions qi{x,s) over x, down to Xmin = ^lrI^i summed 
over all quarks. The total x-nucleon cross-section due to the neutralino gluon 
scattering can be similarly obtained integrating the parton cross-section with 
the gluon structure function of the nucleon. 

The soft breaking terms of the MSSM are characterized by the following basic 
parameters [jT^: the masses of i scalar fields, mg, at the GUT scale, the masses 
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of j gaugino fields, m^^^y ^Iso at GUT scale, the ratio of the two vev's of the 
Higgs fields tan/3 = V2/V1, and the Higgs mixing parameter /i. To proceed, 
we choose two different scenarios. In the first one, we assume universality of 
scalar masses mo and fermion masses mi/2 the GUT scale. In this case 
the neutralino is gaugino dominated in most part of the parameter space 
of the MSSM. To obtain a lower bound for the cross-sections we use the 
configuration with the largest values of parameters compatible with a no- 
fine tuning condition fl^ (mo = 308 GeV, mi/2 = 390 GeV, and /i = 561). 
We fix tan/3 = 8 throughout. This corresponds to a neutralino with mass 
= 160 GeV and a gaugino part of P = Zf^ + = 0.99. The masses 
of the squarks are M^ji ^ 1000 GeV, except the lightest stop which has 
M ^ 780 GeV. 

In the second scenario, we break universality for the two Higgs doublets, 
keeping the universal value mo for all other scalars and mi/2 for gauginos. 



As a result mixed and Higgsino dominated configurations appear [0. We 
choose the configuration with minimal gaugino part, which is obtained for 
mo = 1204 GeV, mi/2 = 295 GeV, /i = 108. The masses of the squarks of 
the first two generations are Ml,r ~ 1400 GeV, while masses of the third 
generation range from 1400 — 730 GeV. 

The resulting cross-sections are shown in Fig. 3 for the case of universality 
(gaugino dominated neutralino), and in Fig. 4 for a non-universal case and 
a higgsino-dominated neutralino. In each figure, the cross-section of gluon- 
neutralino scattering (solid line) and the cross-sections of resonant quark- 
neutralino scattering (dashed lines) for Ftot = zT^qL^u ^ Q + x) with z = 1 
(top) and z = 10 (bottom) are shown. The cross-sections of all subprocesses 
start to grow at energies s ^ 10^ (GeV)^ ~ r- '^^^ ^^^^ with s is caused 
by the decrease of Xmm = Ml j^/s and = (^l,_r + ^qY/s, and the corre- 
sponding decrease of the number of partons with sufficient momentum in the 
nucleon. If squarks do not decay mainly into neutralino, i.e. Ftot ^ ^{(1l,r. — > 
q + x), neutralino-gluon scattering gives in both cases the dominant contribu- 
tion to the total cross-section. At energies s ~ 10^° (GeV)^ or ii^^ ^ 5 ■ lO^'^'^eV, 
the neutralino-nucleon cross-section is about 10^^^ — 10~^^ cm^, i.e. slightly 
lower than the neutrino-nucleon cross-section. 

Let us consider now x + e — > e — > hadrons which is similar to the Glashow 
resonant scattering t'e+e W ^ /i+^'/x- The resonant energy of the neutralino 
is 

= M|/(2me) = 9.8 ■ 10^(Me/10^ GeV)^ GeV (5) 



and the cross-section is also given by Eq. (|). Now, Ftot is the total decay 
width of e determined by the decay channels to electron and neutralino and 
neutrino and chargino and M is the selectron mass. For the frequency of events 
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produced by the neutralino flux IxiE) in a detector with Ne target electrons 
one easily obtains 

u = nN,EoIx{Eo)aes, (6) 

where Q is the solid angle and 

(Jeff = ~ 2 • 10"^^ cm^ (7) 

is the effective cross-section. The factor k depends on the composition of the 
neutralino and is assumed to be of order one. The resonant events with fre- 
quency (H) are produced as a narrow peak at and give an unique signature 
for neutralino. 



5 Interactions of (7-hadrons and its status as LSP 

The interaction of UHE ^-hadrons (G) was already considered in some detail 
in Ref. [|^ for the case of glueballino. Two values determine the interaction of 
a UHE ^-hadron with a nucleon. The flrst one is the radius of the ^-hadron. 
This radius is inversely proportional to the reduced mass of the system and 
estimates of Ref. ^ give for the total GN-cross-section cr ~ 1 mb. In Ref. 0, 
it is argued that this cross-section is of order Aq^^ ~ 10 mb. However, for pro- 
duction of EAS in the atmosphere only interactions with large energy transfer 
are effective. For a gluino with mass rrig ~ 3 GeV or rrig > 150 GeV the large 
energy transfer corresponds to scattering with large and thus to small 
cross-section. 

Now we consider the diffractive interaction of UHE ^-hadron with nucleons. 
Then ^-hadrons exchange with nucleons the 4-momentum Q in the t-channel 
fragmentating into jets of hadrons (including ^-hadron). For a given energy 
transfer y = (E - E')/E, 

QLn = M'y{{l-yr'-myM'), (8) 

where M > itlq is the invariant mass of the fragmentation jet. Unless y is very 
small, Qmin is large and the corresponding cross-section a{y) ~ l/(5^jn(?/) is 
small. We conclude thus that a heavy ^-hadron (m^ ~ 3 GeV or tjiq > 
150 GeV) behaves in the atmosphere like a penetrating particle, while a very 
light ^-hadron interacts like a nucleon. 

Let us discuss now the status of the gluino as LSP. Accelerator experiments 
give the lower limit on the gluino mass as rrig > 150 GeV 0. The upper limit of 
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the gluino mass is given by cosmological and astrophysical constraints, as was 
recently discussed in 0. In this work it was shown that if the gluino provides 
the dark matter observed in our galaxy, the signal from gluino annihilation 
and the abundance of anomalous heavy nuclei is too high. Since we are not 
interested in the case when gluino is DM particle, we can use these arguments 
to obtain an upper limit for the gluino mass. Calculating the relic density of 
gluinos (similar as in Q) and using the condition Qg <C f^cDM, we obtained 
< 9 TeV. 

Now we come to a very strong argument against the existence of a light stable 



or quasistable gluino [|T8|. It is plausible that the glueballino (gg) is the lightest 
hadronic state of gluino ^,0. However, gluebarino, i.e. the bound state of 
gluino and three quarks, is almost stable because baryon number is extremely 



weakly violated. In Ref. |0 it is argued that the hghtest gluebarino is the 
neutral state (guud). These charged gluebarinos are produced by cosmic rays in 
the earth atmosphere |TB|, and light gluino as LSP is excluded by the search for 



heavy hydrogen or by proton decay experiments (in case of quasistable gluino). 
In the case that the lightest gluebarino is neutral, see [|I| , the arguments of 
still work if a neutral gluebarino forms a bound state with the nuclei. Thus, a 
light gluino is disfavored. 

The situation is different if the gluino is heavy, rrig > 150 GeV. This gluino can 
be unstable due to weak R-parity violation |^ and have a lifetime Tg > 1 yr, 
i.e. long enough to be UHE carrier from remote parts of the Universe. Then 
the calculated relic density at the time of decay is not in conflict with the 
cascade nucleosynthesis and all cosmologically produced ^-hadrons decayed 
up to the present time. Moreover, the production of these gluinos by cosmic 
rays in the atmosphere is ineffective because of their large mass. 



6 Discussion and Conclusions 



We estimated the fluxes of UHE LSP produced by the decays of supermassive 
X-particles {nix > 10^'^ GeV). These particles can be Big Bang relics or 
can be produced by topological defects. We performed a simple Monte-Carlo 
simulation for the development of such a cascade. The fluxes of UHE LSP are 
numerically evaluated for two cases: cosmic necklaces as example of topological 
defects, and supermassive Big Bang relic particles. 

Most attention is given to neutralino, as LSP. We calculated the x-^-cross- 
sections at very high energies for two versions of "standard" MSSM with soft 
breaking terms: with universal scalar mass term and non-universal one. The 
typical values of cross-sections at extremely high energies are a ~ 10"'^^ cm^ 
and thus the detection of UHE neutralinos is a difficult task, which needs a 
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special discussion. If the masses of the squarks are near their experimental 
bound Ml^r ~ 180 GeV, the cross-section can be 60 times higher. 

What is important to realize is that neutralinos, which are probably the most 
abundant form of matter in the Universe, can naturally exist in the form of 
high-energy particles. The signature of UHE neutralino is a narrow resonance 
peak at energy E ~ 10^ GeV(Me/l TeV)^ (see Eq.(13)). 



As another example of LSP we considered the gluino. This is a hypothetical 
case, because in all elaborated SUSY models the gluino is not the LSP. Apart 
from a controversial mass window, 1 — 4 GeV, the gluino mass is limited from 
below by accelerator experiments as Trig > 150 GeV, while from above it is 
limited cosmologically [rrig <C 10 TeV). The light gluino as LSP is further 
disfavored due to the upper limit on concentration of heavy hydrogen or by 
searches for proton decay [|l^. The heavy unstable gluino with mass rrig in 
the interval 150 — 1000 GeV and the lifetime > 1 yr is free from the above 
constraints and can be the carrier of UHE signal from remote parts of the 
Universe. Heavy ^-hadrons behave in the earth atmosphere like penetrating 
particles: they have large cross-section, cr ~ 1 mb, for very small energy trans- 
fers but very small cross-sections for large fractions of energy transfer y > 0.1. 
The light ^-hadrons interact like the usual nucleons. 



The energy spectrum of light ^-hadrons has the usual GZK-cutoff, while for 
heavy ones the energy of the GZK-cutoff is shifted to larger energies as shown 
in Fig. 2. In conclusion, we think that gluino is disfavored both as LSP and 
the carrier of UHE signal. 



Acknowledgement 



We acknowledge with gratitude the participation of A. Vilenkin in the initial 
stage of this work. We are grateful to V. De Alfaro, A. Kaidalov and V. 
Matveev for useful discussions and to M. Nagano for providing us with the 
experimental data used in Fig. 1 and 2 and for helpful remarks. M.K. was 
supported by a Feodor-Lynen scholarship of the Alexander von Humboldt- 
Stiftung. 



Note added: At the moment of the submission of this paper, a paper by J. 
Adams et al. ( |hep-ex / 9 7090 28| ) has appeared which excludes a light glueballino 
in the mass and life time ranges of 1.2 — 4.6 GeV and 2 • 10~^° — 7 ■ 10~^ s. 
While these data are not consistent with the model of Farrar et al. [|l],3, they 
do not exclude a long-lived gluino with r > 1 yr which we critically discussed. 
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FIG. 1: Predicted fluxes from decaying X-particles with rx = Q ■ 10~^^: neu- 
trinos (curve neutralinos (curve I^) and g-hadrons from the halo and the 
extragalactic space. The data points are based on the compilation made in Ref. 
[16]. 




FIG. 2: Predicted fluxes from cosmic necklaces with r^//, = 2 • lO^'^ GeV^: 
neutrinos (curve /y), neutralinos (curve and g-hadrons with mass wIq = 
2.5 GeV, niQ = 2.0 GeV and = 1.5 GeV. 
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FIG. 3: Neutralino-nuclcon cross-sections as functions of s = 2E^m]\T for a 
gaugino dominated neutralino. The dashed hnes correspond to the resonant 
subprocess x + q — )■ qL,R all with Ftot = zr{q — )■ q + x) and z = 1 and 
10, respectively. The solid curve represents the subprocesses X + 9 ^ + 
iLA^L,R)- 
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FIG. 4: Ncutralino-nuclcon cross-sections as functions of s = 2E^m]\T for a 
higgsino dominated neutralino. The dashed lines correspond to the resonant 
subprocess x + q — )■ qL,R all with Ftot = zr{q — )■ q + x) and z = 1 and 
10, respectively. The solid curve represents the subprocesses X + 9 ^ + 
iLA^L,R)- 



